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1.  Introduction. 


Hotelling  (1939)  and  Weyl  (1939)  initiated  a  profound  line  of  geometric  research  (e.g. 
Griffiths,  1978,  Langevin  and  Shifrin,  1982,  Gray,  1982),  motivated  by  the  following  statistical 
question.  Suppose  y,  =  /3 /,(0)+£t  (i  =  1,2 ,...,«),  where  the  /,  are  known  functions  depending 
on  an  unknown,  perhaps  multidimensional  parameter  0,  and  the  e,  are  independent  A^O,  a2) 
errors.  What  is  the  significance  level  of  the  likelihood  ratio  test  of  Ho  :  0  =  0?  (One  can 
also  consider  the  more  general  model  y;  =  (o,  z,)  +  /3 /,(0)  +  £,,  where  a  is  a  p-dimensional 
parameter  and  the  x,  are  known  p-dimensional  vectors,  as  we  show  in  Section  4.) 

One  of  Hotelling’s  principal  motivating  examples  is  fi(0)  =  cos (pt*  +  w),  where  the  t,  are 
known  constants  and  0  =  (p,u;).  Another  example  is  the  broken  line  regression  /,(0)  = 
( ti  -  0)+ .  Davies  (1987)  gives  an  interesting  discussion  of  each  of  these  problems  from  a 
different  viewpoint. 

It  is  easy  to  see  that  the  likelihood  ratio  statistic  for  testing  H0  :  (3  =  0  is  equivalent  to 

maxIfE/.^jV^WSy,2]]. 

Setting  f(0)  =  ( f\(0 ), . . . ,  fn(0))  and  y  =  (yi, . . .,  yn),  can  write  this  statistic  in  the  form 

raax{(/(«),v)70l/(9)||2»9||s]}. 

In  terms  of  the  unit  vectors  j(0)  =  f{0)/\\f{0)\\  and  U  =  y/||y||,  the  rejection  region  is  of  the 
form 

max{7(0),  U)2  >  u?2. 

6 

The  rejection  region  can  be  described  geometrically  as  the  union  of  two  tubes  in  the  unit  sphere 
5n_1,  one  about  j(0)  and  one  about  —  7(0),  of  geodesic  radius  cos-1(u>).  Here  the  tube  about 
7(0)  of  geodesic  radius  ip  is  the  set  of  all  u  €  Sn-1  within  geodesic  distance  of  the  manifold 
7(0)  as  0  ranges  over  some  parameter  space.  Under  Ho  U  is  distributed  uniformly  on  Sn~l 
and  hence  the  significance  level  of  the  likelihood  ratio  test  is  the  normalized  surface  area  on 
Sn~ 1  of  the  union  of  the  two  tubes. 

The  purpose  of  this  paper  is  to  review  the  Hotelling-Weyl  approach  and  explore  its 
statistical  implications  by  a  detailed  discussion  of  the  two  concrete  examples  mentioned  above. 
Additional  applications  are  described  by  Naiman  (1986a,b),  Knowles  (1986).  and  Johansen 
and  Johnstone  (1988). 

The  theoretical  basis  for  our  discussion  of  the  broken  line  regression,  which  involves  a  one 
dimensional  parameter  0,  is  given  by  Naiman  (1986)  and  Johnstone  and  Siegmund  (1987).  The 
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theoretical  situation  for  multidimensional  5  is  less  satisfactory.  Since  Weyl’s  theorem  concerns 
manifolds  without  boundary  whereas  many  statistical  problems  and  Hotelling’s  example  in 
particular  involve  manifolds  with  boundary,  we  begin  by  giving  a  version  of  Weyl’s  theorem  for 
two  dimensional  manifolds  (surfaces)  with  boundary  embedded  in  5”"1.  The  two  dimensional 
case  is  especially  interesting  because  a  term  which  would  be  onerous  to  evaluate  numerically 
can  be  combined  with  other  terms  and  evaluated  analytically  by  an  application  of  the  Gauss- 
Bonnet  Theorem.  Naiman  (1987)  discusses  a  special  class  of  d-dimensional  manifolds  with 
boundary,  but  even  in  the  case  d  =  2  his  result  is  somewhat  differently  formulated. 

Our  numerical  studies  have  two  purposes,  (i)  Since  the  Hotelling-Weyl  theory  is  exact 
only  for  tubes  of  small  radii,  i.e.  for  small  significance  levels,  we  indicate  through  these  concrete 
examples  the  accuracy  provided  by  the  theoretical  results  at  conventional  significance  levels, 
(ii)  Since  the  theoretical  results  contain  coefficients  depending  on  the  sample  size  and  dimension 
of  the  parameter  space,  which  can  be  difficult  to  compute,  we  investigate  numerically  the  effect 
of  approximating  or  neglecting  some  of  these  coefficients. 

The  paper  is  organized  as  follows.  Section  2  contains  our  version  of  Weyl’s  theorem  for 
surfaces  with  boundary.  Once  the  theorem  is  formulated  the  proof  is  not  difficult  for  the 
careful  reader  of  Weyl’s  paper.  Nevertheless,  we  give  the  argument  in  some  detail  because 
the  restriction  to  surfaces  makes  possible  a  completely  elementary  (in  the  sense  of  differential 
geometry)  exposition,  which  we  hope  will  make  the  subject  more  widely  known  by  statisticians. 
Section  3  specializes  our  general  theorem  to  the  problem  of  testing  for  a  periodic  term  and 
compares  some  numerical  results  with  the  outcomes  of  a  simulation  experiment.  Section  4 
discusses  the  problem  of  testing  for  a  break  in  a  regression  line.  For  the  most  part  Sections  3 
and  4  can  be  read  independently  of  Section  2.  In  Section  5  we  try  to  draw  some  qualitative 
conclusions  and  discuss  possible  extensions  of  our  results. 

2.  The  Volume  of  Tubes  about  Surfaces  with  Boundary. 

Suitable  references  for  the  following  developments  are  Millman  and  Parker  (1977)  and  Do 
Carmo  (1976).  The  former  is  notationally  consistent  with  Weyl  (1939)  and  with  modern  usage, 
and  consequently  that  notation  is  used  here.  The  latter  gives  a  more  thorough  and  often  more 
adequate  treatment  of  the  concepts  needed  below. 

As  in  Hotelling  (1939)  and  Weyl  (1939)  we  consider  first  the  technically  simpler  case  of 
surfaces  M  in  IRn  and  then  indicate  the  modifications  required  for  surfaces  in  the  unit  sphere 
5n_1  C  mn. 

We  assume  that  a  regular,  oriented  surface  M  C  Hln  is  given  locally  by  a  C3  function 
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x  =  i(u1,«2)  defined  on  an  open  set  in  IR2.  Assume  also  that  the  boundary  of  M,  say  dM , 
is  given  by  a  piecewise  regular,  positively  oriented  curve  7,  parameterized  by  arc  length.  (The 
reader  who  wishes  to  avoid  the  technical  concept  of  orientation  may  consider  the  special  case 
in  which  M  is  given  globally  by  a^u^ti2)  .) 

The  tube  about  M  of  radius  a  is  the  set  of  all  points  whose  minimum  distance  to  M  is 
less  than  a.  We  say  that  overlap  occurs  at  a  point  y  in  the  tube  if  there  exists  more  than  one 
point  in  the  manifold  which  is  closest  to  y. 

Let  \M |  and  \dM\  denote  the  surface  area  of  Af  and  length  of  dM  respectively.  Let 
denote  the  volume  of  the  unit  ball  in  lRn  and  u)n  the  volume  (surface  area)  of  the  unit  sphere 
Sn  ClRn+1  (fl„  =  jrn/2/r(n/2  +  l),wn_i  =  2W2/r(n/2);  our  is  Weyl’s  w„). 

Theorem  1.  Assume  that  the  exterior  angles  at  the  vertices  (if  any)  of  dM  are  positive  in 
the  sense  that  the  tangent  to  7  rotates  through  a  positive  angle  at  each  vertex.  For  all  a 
sufficiently  small  that  no  overlap  occurs  in  the  tube  of  radius  a  about  A/,  the  volume  V(a)  of 
the  tube  is  given  by 

V(a)  =  |M|fln_2an"2  +  2-1|0Af|fin-1a"-1  (1) 

+  2irx(M)ft„-2an/n> 

where  x(M)  is  the  Euler-Poincare  characteristic  of  M  . 

Remarks,  (i)  the  factor  2-k\{M)  in  (1)  arises  indirectly  as 

/  KdA  +  f  kgds  +  E0,,  (2) 

Jm  JdM 

where  K  is  the  Gaussian  curvature  and  dA  the  element  of  surface  area  of  A/,  kg  is  the  geodesic 
curvature  and  ds  the  element  of  arc  length  of  dM ,  and  the  9,  are  the  angles  of  rotation  of  the 
tangent  to  7 (s)  at  the  vertices  of  dM.  That  (2)  equals  2k\{M)  is  the  Gauss-Bonnet  Theorem 
(Do  Carmo,  1976,  p.  274).  Since  numerical  evaluation  of  the  individual  terms  in  (2)  would  be 
onerous  whereas  x{M )  is  easily  determined  and  in  statistical  problems  is  often  zero  or  one,  this 
simplification  of  (2)  is  quite  fortunate.  That  the  various  ingredients  of  (2)  are  involved  in  Y(a) 
is  not  surprising.  That  they  all  have  the  same  coefficient  and  hence  can  be  combined  appears 
in  our  proof  as  an  accident  of  calculation  for  which  we  do  not  have  a  satisfactory  geometric 
explanation. 

(ii )  The  condition  that  no  overlap  occurs  is  both  a  local  and  a  global  condition.  The  tube  over¬ 
laps  itself  locally  whenever  one  of  the  Jacobian  matrices  appearing  in  the  following  argument 


has  a  vanishing  determinant.  In  that  case  the  curvature  of  the  surface  is  so  great  that  the 
mapping  which  parameterizes  points  in  the  tube  is  not  locally  one-to-one.  The  occurrence  of 
global  overlap  is  not  so  easily  characterized.  If  no  local  overlap  occurs,  the  right  hand  side  of 
(1)  is  an  upper  bound  for  the  left  hand  side.  Johansen  and  Johnstone  (1988)  give  an  algorithm 
to  determine  the  smallest  radius  at  which  overlap  occurs  in  a  tube  about  a  curve. 

(iii)  If  some  of  the  exterior  angles  at  vertices  of  dM  are  negative  and  no  self  overlap  occurs 
except  the  local  self  overlap  at  those  vertices,  V(a )  is  <  the  right  hand  side  of  (1). 

Proof.  We  give  the  proof  for  n  =  4,  which  simplifies  the  necessary  computations  and  except 
for  one  technical  identity  given  at  the  end  of  the  proof  has  all  the  ingredients  of  the  general 
case.  The  method  involves  evaluation  of  the  differential  of  volume  at  a  point  y  in  the  tube, 
say  dV(y)f  followed  by  integration  over  the  tube.  We  consider  separately  the  case  of  points  y 
which  are  nearest  an  interior  point  of  M  and  points  nearest  to  a  boundary  point.  Initially  we 
follow  Weyl  closely. 

For  points  y  nearest  to  an  interior  point  x  =  ^(u1,^2)  we  have  y  =  x  +  £i”(l)  +  f2u( 2), 
where  n(l)  and  n( 2)  are  mutually  orthogonal  unit  normals  to  the  tangent  space  to  M  at 
x(ul,u2).  Then  y  =  y(u\  u2,6,{2)  and 

dV(y)=  \\yi,y2,y3,yi\\duidu'id^d^2,  (3) 

where  y,  denotes  the  partial  derivative  of  y  with  respect  to  the  argument,  and  the  double 
bars  denote  the  absolute  value  of  the  determinant  of  the  Jacobian  matrix  of  the  enclosed 
(column)  vectors.  We  have 

Vi  =  Xi  +  frMl)  +  &M2)  O' =1,2),  (4) 

1/3  =  n(l),y4  =  n(2), 

where  x,  =  dx/du'  and  «,-(*/)  =  dn(i‘)jdu'.  We  can  express  n,(v)  as  a  linear  combination  of 
Xi,X2,n(l),  and  n(2)  by  (Weingarten  equations) 

n,(t/)  = -EjZ.f(i/)xj  +  ...,  (5) 

where  -L\(u)  is  for  our  purposes  defined  to  be  the  coefficient  of  Xj  in  the  indicated  expansion, 
and  +  . . .  indicates  components  orthogonal  to  the  tangent  space  spanned  by  xj  and  x2-  Hence 
writing 


L\(v)  L\(v) 


we  have  in  matrix  form  by  (4)  and  (5) 


(yi,y2)  =  (xi,*2)(j  -  fiL(i)  -  £2X(  2))  + 


so  (3)  becomes 

dV(y)  =  ||(xi,*2)(/  -  fi£(l)  -  ^L(2)),n(l),n{2)\\du1du2d^2  (6) 

=  lki,*2,n(l),n(2)||  ||/  - &L(  1)  -  £2Z(2)||du1du2£1d62. 


From  the  determinantal  identity 

|«i,  -  - v*|2  = 


we  see  that  the  first  factor  on  the  right  hand  side  of  (6)  equals  g 1^2,  where  (gij)  =  ((x i,Xj)) 
is  the  matrix  of  the  first  fundamental  form  of  M  and  g  =  \gij\.  Since  dA  —  gx!2dv}du2  is  the 
element  of  surface  area  on  M,  we  see  from  (6)  that  the  volume  associated  with  points  y  closest 
to  interior  points  of  M  is 

Vi(a)  =  I  [[  \I-^L(l)-^L{2)\d^2dA.  (7) 

Jm  <** 


Expansion  of  the  determinant  in  (7)  yields 

1  -  fctrl(l)  -  &trl(2)  + |I(1)|  +  Cll^(2)l 
+  66{Xj(l)^(2)  +  Xl(2)L22(l)  -  2L\(l)L\(2)}, 


so 


Vk( 


.)-/  / 
J  M  Jt 


(l  +  ^|£(l)|  +  ^|X(2)|K,d^ 


=  7Td2 


|M|+/  (|I(1)|  +  |I(2)|)/  tfdhdhdA 

=  7ra2|M|  +  4-17ra4  /  (|L(1)|  +  |L(2)|)<M. 

Jm 


To  identify  |Z(1)|  +  |Z,(2)|  as  the  Gaussian  curvature  A',  we  put  artJ  =  dx/du'duJ  and  in  the 
customary  way  (e.g.  Do  Carmo,  1976,  p.  232)  define  Christoffel  symbols  F*  and  coefficients 
Lij{v)  by 

x,j  =  EfcffjXfc  +  'Z,vLlJ{i/)n(v). 


Since  ( ))  =  0,  Xu(i')  =  (xtJ,  n(i/)}  =  -{xt,  n;(i/)),  and  hence  by  (5)  LXJ(v)  =  Ei 
One  can  now  proceed  line  by  line  as  in  standard  proofs  of  Gauss's  Theorema  Egregium  (e.g. 
Do  Carmo,  1976,  pp.  233-34;  Millman  and  Parker.  1977,  pp.  142-43)  to  show  that  the 
intrinsic  evaluation  of  K  in  terms  of  the  Christoffel  symbols  and  their  partial  derivatives 


equals  |I(1)|  +  |I(2)|.  Hence 


Vi(a)  =  3ra2|M|  +  4  1zai  f  KdA.  (8) 

JM 

Now  assume  that  the  closest  point  in  M  U  dM  to  y  is  a  point  7 (s)  €  dM  where  7 (s)  = 
dy/ds  exists.  Then 

y  =  7  +  r/n(0)  +  &n(l)  +  ^«(2),  (9) 

where  n(0)  is  a  unit  vector  in  the  tangent  plane  of  M,  but  orthogonal  to  7  and  as  before  n(l) 
and  n( 2)  are  orthogonal  to  the  tangent  plane.  We  choose  n(0)  to  point  towards  the  interior  of 
the  manifold  (rotated  by  -fx/2  from  7),  so  the  geodesic  curvature  kg  =  kg(s)  =  (7,n(0)),  and 
T)  €  (-a,0]. 

The  volume  element  is 


dV(y)  =  \\yi,y2,y3,y*\\dsdridZid{2, 

where  as  before  y,  is  the  partial  derivative  of  y  =  y(s,  t/,£i,£2)  with  respect  to  the  ith  argument. 

Putting  q  =  dq/ds,  we  have  the  Frenet-Serret  equations  (cf.  Do  Carmo,  1976,  p.  261) 

7  =  kgn(  0)  +  kin(l)  +  fc2n(  2) 
n(0)  =  -kgi  +  nr?(l)  +  r2rc(  2) 
n(l)  =  -k\i  -  u(0)  +  t3ti(2) 

n(2)  =  -fc2 7  -  r2n(0)  -  T3n(l), 

where  except  for  kg  the  values  of  the  coefficients  ultimately  will  not  concern  us.  Hence  by  (9) 

dV  =  ||7  +  77h(0)  +  C1h(l)  +  ^n(2),n(0),n(l),n(2)||dsd7?d^2 
=  ||7(1  -  Tjkg  -  frfci  -  {2fc2),n(0),n(l),n(2)||dsd77^1cf^2. 

The  volume  associated  with  these  points  y  is 

vm-LIIL 


rj<0 
n2+<J+^<a2 


(1  -  r]kg  -  krfi  -  k2^2)d^id^2dvds 


=  /  /  7r(a2  -  T72)(  1  -  kgT])drids 

JdM  J-a 

=  27r|9A/|u3/3  +  4_17ra4  j  kgds. 

JdM 


(10) 


The  contribution  to  V(a)  from  points  in  the  tube  which  are  closest  to  vertices  of  dM  is 
clearly 

V3(a)  =  n4a4L6,/2*  =  4_,jra4S«„  (11) 


where  the  O' s  are  the  angles  of  rotation  of  the  tangent  to  y  at  the  vertices. 

Addition  of  (8),  (10),  and  (11)  and  an  appeal  to  the  Gauss-Bonnet  Theorem  as  indicated 
in  Remark  (i)  following  the  statement  of  the  theorem  complete  the  proof  in  the  special  case 
n  =  4. 

The  proof for  general  n  is  essentially  the  same.  The  coefficients  of  KdA,  JgM  kgds, 
and  E0,  are  respectively 


«”/•••/  tfd6,...d£n-2,-an  /...  / 

J  K’_S<1  J  •/„2+{?+n 


vdrjd^i . .  .d^n--'. 


and  (2?r)  1anfln,  which  some  calculus  shows  are  all  equal  to  n_1a"nn_ 2  .  I 

To  prepare  for  the  spherical  case  suppose  y  =  y(ti* , . . tin_1)  is  a  coordinate  patch  of  a 
hypersurface  inlR”.  Let  p  =  {y,y)  and  z  =  y/p  £  Sn_1.  Assume  z  is  one-to-one. 

Lemma  1  (Weyl):  The  volume  element  at  z  £  Sn~l  is  given  by 


dV(z)  =  \}y,yu---,yn-i\\du' ..  .dun~l/pn. 


where  y,  ==  dy/du'. 


Proof.  Since  z  is  the  unit  normal  to  the  hypersurface  z(u , . . ’), 


where  z,  =  dz/du'.  Since 


dV(z)  =  ||z,z1,...,zn_1||du1  ...dun  \ 


z,  =  p  }y,  +  ydp  1  / <9u\ 


by  colli nearity  we  have  |z,zj, . .  .,zn_i|  =  |y,  yx, . . . ,  y„_i |/pn.  I 

Now  suppose  M  C  5n_1,  the  unit  sphere  ini??”.  The  tube  about  A/  of  geodesic  radius  <p 
is  defined  as  the  set  of  all  z  £  S”-1  such  that  ( z,x )  >  cosys  for  some  x  £  M. 

Theorem  2.  Assume  that  M  C  S”-1  and  the  exterior  angles  at  the  vertices  of  d\I  are 
positive.  For  all  yj  sufficiently  small  that  no  overlap  occurs  in  5n_1  in  the  tube  of  geodesic 
radius  y?  about  A/,  the  volume  V(<p)  of  the  tube  is  given  by 

V(yj)  =  (n  -  3)_1wn_4  [|A/j  cosy?(sin  y:)”-3 

+2*x(A/)  r(sinw)n-2du.-]  +  [2(n  -2)]“1u>n_3|0A/|(sinys)n"2. 

Jo 


Pi!' 


si 
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where  \( M )  is  the  Euler-Poincare  characteristic  of  M . 


Proof.  Suppose  z  is  a  point  in  the  tube  closest  to  an  interior  point  x  of  M .  Put  k  =  n  -  3 
and 

y  =  y(u\u2,(  =  x  +  6n(l)  +  •••  +  tkn(k),  (13) 

where  x  =  ^(u^u2)  €  M  and  the  n(y)  €  S”-1  are  mutually  orthogonal  unit  normals  to  the 
tangent  space  of  M  at  x(vJ ,  u2).  Then  we  have  the  representation 


where 

£1  +  •••  +  £*  <  tan2  <p. 


By  Let  ima  1  and  (13) 


dV(z)  =  ||y,yi,y2,n(l),...,n(^)|| 


dv}du2d£i  . . .  dfjt 
(1  +  f?  +  •••  +  ffc)n/2‘ 


The  proof  of  Theorem  1  together  with  some  calculation  shows  that  the  contribution  to  the 
volume  of  the  tube  from  these  points  is 


+  «)"'2 


—  w„_4  [|A/|  J  cos2w(sinu;),l_4<fu; 

4-  (n  -  3)_1  /  (K  -  l)dA  f  \sinui)n~2cL;] 

Jm  Jo 

=  (n  -  3)-1oJn_4  [|jl/|  cos  <^>(sin  y)n_3  +  /  A'</.4  f  (sin  ^')n~2d^] . 

Jm  Jo 


(14) 


The  reason  for  the  appearance  of  K  -  1  in  an  intermediate  equality  of  (14)  is  that  ( 13)  involves 
k  =  n  —  3  of  the  n  —  2  normals  (in  IRn)  to  the  tangent  space  of  M .  The  last  normal  is  r,  the 
normal  to  S"-1  itself.  Its  contribution  to  the  (intrinsic)  Gaussian  curvature  of  M  equals  the 
Gaussian  curvature  of  Sn~l .  which  is  one.  Hence  from  the  n  —  3  normals  represented  in  (13) 
we  obtain  K  —  1. 


For  points  z  closest  to  a  point  of  dM  where  7(.s)  exists,  we  have 


% 


» 


y  =  7  +  vn(0)  +  (it(i  )+••■  +  &«(*)•  2=  y/(  1  +  rf  +  m  +  . .  •  +  ik  )'1". 

where  n( 0)  is  orthogonal  to  i  but  in  the  tangent  spare  of  M  (pointing  into  M  ).  and  n(  1 ). 


y 

v-. 


T)(k)  are  as  before.  Also  rj  <  0  and  r?2  +  f 2  +  . . .  +  <  tan2  9.  By  Lemma  1 


<*V(z)  =  ||y,yi,  •  • . ,  yn-\\ 


dsdr)d£i . . .  d£k 

(1  +  r/2  +  ti  +  ...  +  tf)nl2 


~  Ill,')'  +  rfn(O)  +  S^n(i),n(0),n(l),...,n(*)|| 


dsdrjdt  1 . .  .<i& 

(l  +  *72  +  {? +  ■  -  +  ^)n/2‘ 


Using  the  Frenet-Serret  equations  as  in  the  proof  of  Theorem  1,  we  obtain 

„  r,  ,  \  dsdr)dix  ...d£k 

di  (2)  = (1  - 

where  fcg  =  (7,n(0))  =  -(7,n(0))  is  the  geodesic  curvature  of  q,  and  the  exact  values  of  the 
coefficients  c,  need  not  concern  us.  Hence  the  contribution  to  the  volume  of  the  tube  arising 
from  these  points  is 


V2(<p) 


JdM  I 


(1  -  kgr,) 


*72  +  €j4  -f  <t*n2  v 


dr)d£\  . .  .d£kds 
(i  +  »72  +  tf +  •"•  +  £?  )n/2 


=  [2(n  —  2)]  1un^3\dM |(sin^)n  2  +  (n-3)  ’u;n_4  [  kgds  [  (sinu>)n  2cLj. 

JdM  Jo 

Finally,  the  contribution  arising  from  points  nearest  to  a  vertex  of  dM  is 

V3M  =  (n  -  3)~1u;n_4(2^)  r(sinwj)"-2<L>.  (16) 

Jo 

Addition  of  ( 14 ).  (15).  and  ( 16)  and  an  appeal  to  the  Gauss-Bonnet  Theorem  as  in  Remark 
(i)  following  the  statement  of  Theorem  1  complete  the  proof.  I 

Corollary  1.  Under  the  conditions  of  Theorem  2,  if  U  is  uniformly  distributed  on  S"-1, 


/*{  sup  (x(u',  n2),  T)  >  u'}  = 


2  1  r ( n / 2 )  1  A/ 1  u- ( 1  -  u'2 )  2' 


- -  1  »3/>r[(n-n/2] 

+<4»na.ui(i  -  /,(I  - 

tr1/2r  [( n  -  1)/2J  Jw 

Proof.  The  corollary  follows  from  the  fact  that  the  indicated  probability  is  the  normalized 
volume  of  the  tube  of  geodesic  radius  #  -  cos-1  u\  i.e. 

ncos-1  u’)/u,’„_i . 


3.  Testing  for  an  Harmonic. 


Assume  yi  =  /?cos(pt,-  -fu)  +  e,-  (i  =  1,2, where  the  t,  are  known  constants  and 
th''  f:  are  independent  JV(0,cr2).  The  parameters  0,n,u>,  and  a2  are  all  unknown,  and  we  wish 
to  test  Ho  :  0  =  0  against  Hi  :  0  ^  0  .  It  is  convenient  to  write 

0  cos(mU  +d)  =  01  cos(/it,-)  +  02  sin (/**,-), 

where  0\  —  /?cosu>,/?2  =  —  /Jsinw,  and  consider  the  equivalent  problem  of  testing  Ho  '■  0\  - 
02  =  0.  A  calculation  shows  that  the  likelihood  ratio  test  has  a  rejection  region  of  the  form 

sup  [(71(p),f/)2  +  (72(/i),f/)2]1/2  >  tn,  (17) 

where  7 '(n)  €  5n_1,  (71(/i),72(p))  =  0  for  all  p,  and  U  =  J//||y||  is  uniformly  distributed  on 
Sn_1  if  Ho  is  true.  In  the  special  case  that  the  ti  are  equally  spaced,  the  natural  range  of  /i  is 

(0,Jr). 

Johnstone  and  Siegmund  (1987)  exploit  the  one  dimensional  structure  of  (17)  and  obtain 
by  a  non-geometric  argument  an  upper  bound  for  the  signficance  level  of  the  likelihood  ratio 
test.  We  shall  show  as  an  application  of  Corollary  1  that  their  inequality  is  an  equality  when 
the  significance  level  is  sufficiently  small.  We  also  show  by  a  Monte  Carlo  experiment  that 
when  the  t,  are  equally  spaced  these  analytic  results  provide  very  good  approximations  for  the 
true  significance  level  over  a  broad  range  of  values. 

Note  that 

[(71(/x),t/)2  +  (72(/r),f/)2]1/2  =  sup  [(71,f/)cosu;+ (72,t/>sinu>] 

0<u/<2tt 

=  sup  [(71  cosu-' +  72  sinw,  {/)]. 

0<u><2tt 

Hence,  if  we  put  x(p,co)  =  71(p)cosu>  +  7z(p)sinu>,  Corollary  1  gives  us  the  probability  of 
(17)  in  terms  of 


and  x(A/). 

Since  |j<9x/<9w||  =  1  for  all  p,  \dM\  =  4x.  It  is  easy  to  see  that  \{M)  =  0.  For  example, 
one  can  cut  the  surface  along  w  =  0  to  obtain  a  new  surface  which  clearly  has  Euler-Poincare 


characteristic  equal  to  one  and  a  connected  boundary  with  four  90°  exterior  angles.  Comparison 
of  the  Gauss-Bonnet  Theorem  applied  to  the  cut  and  uncut  surfaces  shows  that  y(A/ )  =  0. 


To  obtain  an  expression  for  |A/|,  observe  that  =  dx/dui  =  -71  sinu;  +  72cosu;  and 
x ^  —  dx/dp  =  7  cosw+7  sinu,  so  ||arw||  =  1  and  ||xj|  =  ||7  cosu;+7  sinu;||.  Since  (7\72)  = 
0,  (i1^2)  =  -(71, 72)  and  hence  (x^^x^)  =  (7*,  72)  cos2  u>  -  (7 \  72)  sin2  u>  =  (71 , 72). 

Substitution  of  these  results  into  the  formula  given  in  Corollary  1  yields 


Corollary  2.  If  7'(/i)  €  Sn  1  (i  =  1,2)  are  non-selfintersecting  regular  C3  curves  and 
(71(^),  72(/i))  =  0  for  all  /ij  <  n  <  /i2,  then  for  all  u>  sufficiently  close  to  1, 


p{  sup  [(71(Ai)^>2+(7V),tr)J]1/2>^|  =(l-u-2 


,2\(n-2)/2 


r  r  i"1'  — + ^ 

Remark.  Using  quite  different  methods,  Johnstone  and  Siegmund  (1987)  showed  that  the 
right  hand  side  of  (18)  is  an  upper  bound  for  the  left  hand  side  for  all  0  <  u-  <  1.  In 
the  one  dimensional  case,  i.e.  for  a  tube  about  a  curve,  Johnstone  and  Siegmund  provide  two 
derivations  of  upper  bounds  for  the  tube  volume  which  are  shown  to  be  equalities  when  no  self 
overlap  occurs  in  the  tube.  However,  their  method  does  not  seem  to  yield  the  equality  in  (18). 
and  our  geometric  method  does  not  appear  to  give  their  inequality  for  all  0  <  u-  <  1.  In  fact 
the  problem  of  obtaining  tight  upper  bounds  for  the  volume  of  tubes  of  arbitrary  radii  about 
manifolds  of  dimension  >  2  appears  to  be  a  difficult  problem. 

For  a  numerical  example  suppose  that  the  f,  are  equally  spaced,  so  by  recentering  we  have 
t,  =  i  —  m,  where  m  =  (n  +  1  )/2,  i  =  1, . .  ,,n.  The  model  is  then 

V  =  0\  f\n)  +  -32/2(p)  +  e, 

where  y  =  (yu  . . yn),e  =  (fi . en)J\n)  =  (cosp(l  -  m) . cos p(n  -  m)),  and  /2(/x)  = 

(sin/r(l  -  m),. .  ,,sinp(n  -  m)).  It  is  readily  verified  that  (/' , /2)  =  0,||/1||2  =  n/2  + 
sin(n/x)/(2 sin^i),  |]/2||2  =  n/2  -  sin(np)/(2sinp),  and  hence  the  7’  in  (17)  are  given  by 
7*  =  /7||/*||  (i  -  U2).  Also  (71 , 72)  =  0  =  For  /ij  =  0,p2  =  x  the  integral  in 

(18)  for  our  special  case  reduces  to 


r*/ 2  re2 


/’filial 


!cos2u?  +  ||72||2sin2w 


•  1  -2 

Obtaining  useful  expressions  for  ||7  ||  and  ||7  ||  requires  a  lengthier  calculation.  Some 
details  are  given  by  Davies  (1987),  who  (assuming  o 2  is  known)  discusses  the  likelihood  ratio 
test  from  a  different  viewpoint  and  is  led  to  the  same  integral.  Davies  observes  that  ||7*(/i)||  ~ 
n/(2  •  31/2)  except  for  n  close  to  0  or  x,  which  suggests  use  of  x2n/ 31/2  as  an  approximation 
for  (19).  He  subsequently  modifies  this  suggestion  to  be  more  consistent  with  his  numerical 
calculations.  Although  the  modified  approximation  seems  excellent,  in  the  numerical  example 
given  below  we  have  computed  (19)  numerically. 

For  various  values  of  n  and  tv,  Table  1  gives  the  right  hand  side  of  (18)  and  the  outcome 
of  a  10,000  repetition  Monte  Carlo  experiment  to  estimate  the  probability  on  the  left  hand 
side  of  (18).  The  analytic  result  seems  to  provide  a  very  good  approximation,  although  it  is 
less  satisfactory  for  large  n. 

Table  1 

The  Probability  (18) 


n 

XV 

Right  Hand  Side  of  (18) 

Monte  Carlo  Estimate 

8 

.90 

.092 

.091 

8 

.95 

.018 

.015 

16 

.75 

.073 

.073 

16 

.80 

.022 

.021 

32 

.55 

.186 

.181 

32 

.60 

.058 

.057 

32 

.65 

.014 

.014 

64 

.40 

.355 

.298 

64 

.45 

.082 

.069 

64 

.50 

.014 

.013 

The  term  in  Corollary  1  which  involves  the  Euler-Poincare  characteristic  of  M  is  of  smaller 
order  of  magnitude  than  the  others  as  u;  — *  1.  For  manifolds  of  dimension  >  2  the  analogous 
terms  do  not  in  general  simplify  as  they  do  in  the  two  dimensional  case,  and  they  can  be  quite 
complicated  to  evaluate  numerically.  Hence  it  is  interesting  to  consider  the  accuracy  of  the 
approximation  we  obtain  by  simply  ignoring  those  terms. 

Our  second  example  is  concerned  with  essentially  the  same  manifold  i(/i,to)  =  cosa! 

+  72(//)sinw  as  before,  but  now  to  is  restricted  to  the  range  0  <  to  <  x.  The  area  of  the  new 


manifold  is  one-half  that  of  the  old.  The  length  of  its  boundary  is 

2n  +  2  f  \\il\\dn, 

Jo 

and  its  Euler-Poincare  characteristic  is  one. 

Table  2  gives  a  numerical  example.  The  columns  headed  )M|,|9M|,  and  x{M)  give  the 
values  of  the  respective  terms  in  Corollary  1.  The  final  column  gives  the  results  of  a  10,000 
repetition  Monte  Carlo  experiment. 

Table  2 

Effect  of  Neglecting  the  Gauss-Bonnet  Term 


n 

w 

|M| 

\dM\ 

X(M) 

Total 

Monte  Carlo 

8 

.86 

.025 

.116 

.115 

8 

.94 

.015 

.016 

16 

.72 

.015 

.079 

.082 

16 

.78 

.021 

.024 

32 

.54 

.111 

.028 

.135 

32 

.58 

.057 

.057 

64 

.176 

.043 

.219 

.186 

64 

.45 

.049 

.040 

The  analytic  result  seems  to  be  about  as  accurate  as  in  the  preceding  example.  Neglecting 
the  term  involving  x(M)  would  have  negligible  effect,  but  neglecting  the  term  involving  \dM\ 
would  make  a  substantial  difference.  Since  both  \M\  and  \dM\  may  be  shown  to  be  proportional 
to  n ,  one  might  have  predicted  these  relations  from  analytic  considerations.  It  is  also  interesting 
that  in  this  case  our  theoretical  result  is  not  known  to  be  an  upper  bound;  and  in  three  rows 
the  Monte  Carlo  estimate  actually  exceeds  slightly  the  theoretical  value. 


4.  Testing  for  a  Change-point  in  Linear  Regression. 

Assume  y,  =  0O  +  (JiU  +  -  0)+  +  £,  (i  =  l,...,n),  where  the  t,  €  [0,1]  are  known 

constants  /3o,/?i,  and  (J  are  unknown  parameters,  and  the  e,  are  independent  N( 0,<r2).  The 
problem  of  testing  H0  =  0  =  0  has  been  discussed,  for  example,  by  Hinkley  (1969,  1971), 
Feder  (19 75),  and  Davies  (1987).  Presumably  one-sided  alternatives,  sav  Hi  :  (3  >  0,  are 
scientifically  more  relevant  and  we  consider  only  this  case.  For  our  numerical  examples  we 
take  the  t ,  to  be  equally  spaced:  t,  =  (i  -  l)/(n  -  1),  i  =  1, . . . ,  n. 


We  begin  with  the  special  case  /30  =  Pi  =  0,  which  fits  directly  into  the  framework  of 
Section  1  with  fi(6)  =  ( ti  —  0)+  ( i  —  Since  7(0)  =  /(0)/||/(0||  is  only  piecewise 

smooth,  Hotelling’s  (1939)  results  are  not  directly  applicable.  Naiman  (1986)  and  Johnstone 
and  Siegmund  (1987)  obtain  a  suitable  extension  which  gives  an  upper  bound  for  the  volume 
of  the  tube  of  geodesic  radius  about  7  for  all  0  <  <  tr/2. 

We  consider  a  slight  generalization  of  the  likelihood  ratio  test,  where  we  may  restrict  the 
set  over  which  we  search  for  a  possible  change-point  to  a  subinterval  of  [0, 1],  say  [O,0o].  (See 
James,  James,  and  Siegmund,  1987,  for  the  effect  this  modification  has  on  the  power  of  the 
likelihood  ratio  test  in  a  related  problem.)  The  test  statistic  is  maxo<0<0o  (7(0),  U)  and  the 
length  of  7  is  )7 j  =  j£°  |)7(0)||d0.  According  to  Naiman  (1986)  and  Johnstone  and  Siegmund 
(1987)  for  all  0  <  w  <  1 

p|  max  (7(0),  U)  >  wl  <  — -  -tt  [\l-  x2)^~3^dx  (20) 

\o<e<e0'  h  '  J  ~  iri/2r[(n-  l)/2]  JWK  ’  V  ' 

+(2jt)-,|7|(1-  u'2)("-2)/2. 

Assume  now  that  t,  =  ( i  -  1  )/(n  -  1)  (i  =  l . n).  Some  calculation  shows  that 

||7(0)||  =  0  for  (n-2)/(n-l)  <  0  <  1  and  for  i//(n-l)  <  0  <  (i/  +  l)/(n-l)(i'  =  0, 1, . . .,  n-3) 


where  A„(0)  =  02  +  b^B  +  c„,  bv  =  -(n  +  v)/(n  -  1),  and  c„  =  {2[(n  -  l)2  +  3(tj  -  l)u  +  i/2)]  + 
3(n  +  v  —  1)  +  l}/  [6(n  -  l)2] .  It  is  easily  checked  that  A^(0)  >  0  for  u  <  n  -  3.  so 

r(n— l)-1  [(n-l)0oj  ™  f  („  _  1  _  „}2  _  1  )  r(*'+l)/(n-l) 


||7(0)||d0 


f(n-  1  -  i/)2  —  l\  ' 

ZT'qX  12(n  —  l)2  J  Jv/(n—l) 


dB/A^B ), 


where  m  =  min(n  —  3,  [(n  —  l)0oj )•  Integration  gives 


f(n-i)-’  L(»-i)«oJ 


=  E 


31/2(t,  4.  2  -  Tl) 
{(n  -  1  -  I/)2  -  l}1/2 


—  tan 


31/2(1/  —  Tl)  1 

{(n  -  1  -  u)2  -  1 } 1  /2  J 


With  the  help  of  (21)  one  can  evaluate  the  right  hand  side  of  (20)  numerically.  (Note  that 
the  first  term  is  just  P{t„_!  >  ( n  -  l)1/2tu/(  1  -  te2)1/2},  where  tn  denotes  a  random  variable 
having  Student’s  distribution  with  n  degrees  of  freedom.)  Some  additional  analysis  shows  that 


as  n  — ►  00 


.VaWa'A 


f 

'6 


WOK 


r  ||7(*)||<W  -  -31/2[log(l  -  tf0)]/2,  (22) 
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which  may  be  used  instead  of  (21)  to  evaluate  approximately  the  right  hand  side  of  (20). 

Table  3  compares  the  upper  bound  (20)  with  the  results  of  a  10,000  repetition  Monte 
Carlo  experiment.  Both  (21)  and  (22)  are  used  to  evaluate  the  right  hand  side  of  (20),  and 
both  yield  reasonably  good  approximations  to  the  true  probability. 

Table  3 

Change-point  in  Regression  (0o  =  0\  =  0). 


n 

w 

0o 

Theory  (20) 

Exact(21)  Approx.  (22) 

Monte  Carlo 

5 

.729 

(n  —  2)/(n  — 

1) 

.10 

.111 

.095 

10 

.637 

(n-2)/(n- 

1) 

.05 

.055 

.049 

20 

.476 

(n-2)/(n- 

1) 

.05 

.055 

.047 

20 

.452 

.80 

.05 

.052 

.043 

50 

.274 

(n-2)/(n- 

1) 

.10 

.108 

.090 

50 

.240 

.80 

.10 

.107 

.094 

50 

.317 

(n-2)/(n- 

1) 

.05 

.054 

.048 

50 

.286 

.80 

.05 

.054 

.047 

We  now  drop  the  assumption  00  =  01  =  0.  As  indicated  in  a  slightly  simpler  context 
by  Hotelling  (1939)),  one  can  reduce  this  case  to  the  one  considered  previously  as  follows. 
Writing  our  model  in  the  form  t/;  =  0O  +  0\t{  +  /?/,(#)  +  e,  (i  =  l,...,n)  and  introducing 
the  notation  y  =  (yi  ,...,yn),  t  =  (<i,.  ..,<„),  f(0)  =  (/i(0), e  =  («i, ...,£„),  e  = 
(1, . . .,  1),  y  =  n-1£y,,  t  =  and  f{9)  =  n_1E/;(0),  we  let  y,  =  y,  -  / %  -  0\U>  where 

0i  =  (t  —  <e,y)/||<  -  fe||2  and  0o  =  y  -  0\ t  are  the  maximum  likelihood  estimators  of  0i  and 
0o  under  Hq  :  0  =  0.  It  is  readily  verified  that 

y.  =  0109)  +  e,  (j  =  i,...,n). 


mo)  =  mo)  -  m  ~ 


Hz  _  /^ 


where 


Also  7(0)  =  /(0)/ 1|/(0)||  defines  a  curve  on  the  n  -  3  dimensional  sphere 

,un)  :  Eu-  =  l,£u,  =  0, E(t,  -  t)u,  =  0}, 

and  under  Ho  :  0  =  0,y/||y||  =  U  is  uniformly  distributed  on  this  sphere.  It  follows  that 
(/(0),y)  =  (/(0),ji)  and  (20)  gives  an  upper  bound  for  the  significance  level, 


of  the  likelihood  ratio  test  provided  (7!  is  replaced  by  |7|,  the  length  of  the  curve  7,  and  n  is 
replaced  by  (n  —  2)  to  account  for  the  degrees  of  freedom  lost  in  estimating  0O  and  /?,. 

It  is  now  in  principle  straightforward  to  evaluate  the  arc  length  of  j(0)  and  hence  obtain 
an  upper  bound  for  (23).  There  does  not  seem  to  be  an  exact  expression  as  simple  as  (21); 
and  since  (22)  seemed  to  provide  a  reasonably  good  approximation  in  the  previous  case,  we 
use  the  analogous  approximation  here.  After  substantial  calculation,  one  finds  that  as  n  — <■  00 

r  ||7(0)||d0  -  (31/2/2)log  [0,(1  -  0o)/0o(l  -  0i )] -  (24) 

Je  o 

Table  4  compares  the  analytic  upper  bound  for  (23)  with  the  results  of  a  10,000  repetition 
Monte  Carlo  experiment.  The  analytic  upper  bound  is  given  by  (20)  with  n  replaced  by  n  -  2 
and  (7 (  replaced  by  |7|.  We  have  used  the  right  hand  side  of  (24)  with  0o  =  l/(n  -  1)  and 
=  (n  —  2)/(n  —  1)  as  an  approximation  to  |7j.  For  each  value  of  n  there  is  one  Monte 
Carlo  experiment,  so  the  entries  for  fixed  n  and  different  values  of  tr  are  dependent.  The 
results  indicate  that  except  for  small  n  the  theoretical  approximation  is  reasonable,  albeit 
not  excellent.  Table  3  suggests  that  we  could  perhaps  halve  the  error  of  the  approximation  by 
computing  I7I  numerically. 
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Table  4 

Change-point  in  Regression  (/3o,^i  unknown) 


n 

W 

Theory 

Monte  Carlo 

5 

.963 

.10 

.059 

5 

.989 

.05 

.018 

10 

.693 

.10 

.080 

10 

.763 

.05 

.038 

10 

.867 

.01 

.008 

20 

.495 

.10 

.085 

20 

.553 

.05 

.043 

20 

.656 

.01 

.010 

50 

.321 

.10 

.077 

50 

.359 

.05 

.041 

100 

.232 

.10 

.079 

100 

.259 

.05 

.041 

5.  Discussion. 


In  this  section  we  consider  briefly  some  possible  extensions  of  Theorem  2.  The  most 
obvious  concerns  tubes  about  higher  dimensional  manifolds  (with  boundary).  Unfortunately 
there  does  not  seem  to  be  any  way  in  general  to  avoid  some  nasty  looking  computational 
problems. 

For  example,  suppose  A/  is  a  three  dimensional  submanifold  of  Sn-1  with  a  two  dimen¬ 
sional  boundary  DM.  It  is  possible  for  dM  itself  to  have  a  (one-dimensional)  boundary,  but 
for  simplicity  we  assume  it  does  not. 


The  method  of  proof  of  Theorems  1  and  2  shows  that  the  contribution  to  the  volume  of 
the  tube  of  geodesic  radius  about  A/  coming  from  points  in  .S'"-1  which  are  closest  to  an 
interior  point  of  A/  is 


V,( 


HV7 


i\l 


+i\< tanJ 


l  +  tf(J-3)] 


. 


1  +  +  •  • 


+  Sl)n/2 


dA. 


Here  k  =  n  -  4,d.4  is  the  element  of  surface  area  on  A/,  and  J  is  a  curvature  on  M  defined  as 

follows.  If  (L\(u))(v  =  1 . k+  1 )  are  3  X  3  matrices  defined  as  in  (5)  relative  to  the  normals 

n( v),v  =  1 _ _ k,  to  A/  in  the  tangent  space  of  5n_I  and  the  normal  u(k  +  1 )  to  M  which  is 


also  normal  to  S"  1 ,  then  J  is  the  sum  over  v  —  1, . . k  +  1  of  the  sum  of  the  three  pairwise 
products  of  the  eigenvalues  of  (L^i')),  i.e.  the  second  symmetric  function  of  ( L\(v )). 


V2(<p) 


The  contribution  to  the  tube  volume  from  points  closest  to  dM  can  be  shown  to  equal 

dr)d(\ , . .  .,dj\ 


IdM  /j+<? 


^<0 

+  +<2<I4„j  „ 


[1  -2r,Ha(0)  +  ti(K9-  1)] 


(l  +  r?2  +  tf  +  ...  +  ^)n/2 


dAa, 


where  dAd  is  the  elemenet  of  surface  area  on  dM,  A'a  is  the  Gaussian  curvature  of  dM ,  and 
Hd( 0)  is  the  part  of  the  mean  curvature  of  dM  associated  with  that  normal  n(0)  to  dM  which 
lies  in  the  tangent  space  of  M  (and  points  into  M). 


Addition  of  Vj  and  V2  and  some  calculation  give  an  expression  for  the  volume  of  the  tube, 
which  involves 

/  dA,  f  dAd,  I  JdA  +  2  [  Hd(0)dAd, 

JM  JdM  JM  JdM 

and  fdM  h'sdAd.  The  last  integral  can  be  simplified  by  an  application  of  the  Gauss-Bonnet 
Theorem.  The  others  must  be  evaluated  numerically  in  most  cases.  Numerical  computation 
of  the  third  expression  appears  to  be  rather  complicated  in  general,  although  in  special  cases 
the  integrands  can  be  simplified  analytically. 


We  are  primarily  interested  in  small  probabilities,  and  it  is  easy  to  see  that  here  as 
in  Theorems  1  and  2  the  terms  involving  curvatures  of  M  and  dM  are  of  smaller  order  of 
magnitude  as  w  — *  1  than  those  involving  the  surface  areas,  |A/|  and  \dM\.  Thus  we  might  use 
only  these  comparatively  easily  computed  terms  as  an  approximation.  The  second  example  in 
Section  3  shows  that  such  an  approximation  would  be  quite  good  in  that  case. 

We  record  here  the  first  two  terms  of  the  small  tube  probability  for  an  m  dimensional 
manifold  with  boundary  embedded  in  Sn_1.  We  assume  M  C  Sn_1  is  defined  locally  by 
x  =  x(t),  t  €  T  C  IRm ,  let  |M|  and  \dM\  denote  the  surface  area  of  M  and  dM ,  and  let  U  be 
uniformly  distributed  on  5n_1.  Then  as  tr  —  1 


P{  max(x(/),  V)  >  tr}  =  -  tt,m- 1(1  _  u,2)(n-m-l)/2 

1  I  '  1  >  7idm  +  1)/2r[(n  -  m  +  l)/2] 


(25) 


4~'T(n/2)\dM\ 


jrm/2r[(n  -  m  +  1  )/2] 


7yU'm-2(l  -  U'2)(n"m)/2  +  0((1  -  u-2  )(n-Af  +  l  )/2) 


Complete  development  of  these  ideas  is  a  project  for  future  research.  See  Naiman  (1987) 
for  a  theoretical  beginning  in  this  direction  and  Johansen  and  Johnstone  (1988)  for  a  thought 
provoking  application. 


Naiman  (1986)  and  Johnstone  and  Sicgmund  (  1987)  show  that  Hotelling's  formula  giving 
the  exact  volume  of  a  tube  of  small  radius  about  a  non-closed  curve  actually  provides  an 


upper  bound  for  the  volume  of  a  tube  of  arbitrary  radius.  (The  result  also  applies  to  closed 
curves  provided  they  are  given  fictitious  end  points  and  regarded  as  not  closed.)  It  would  be 
interesting  to  obtain  good  upper  bounds  for  the  volume  of  tubes  of  arbitrary  radii  about  higher 
dimensional  manifolds.  However,  simple  examples  of  surfaces  in  IR 3  show  that  the  situation  is 
much  more  complicated  than  for  curves.  In  particular,  for  a  cylinder  of  unit  radius  and  height 
h  the  formula  for  the  volume  of  tube  of  radius  a  <  1,  to  wit  Airah  +  2?r 2a2,  does  not  give  an 
upper  bound  for  the  volume  of  a  tube  of  radius  a  >  1  when  h  is  large.  On  the  other  hand  the 
formula  for  the  volume  of  a  tube  of  radius  a  <  1  about  S2  C  IR 3,  87ra(l  +  a2/ 3),  provides  an 
upper  bound  for  the  volume  for  all  a  although  for  S1  C  IR 2  the  closed  curve  must  be  given 
fictitious  endpoints  in  order  to  obtain  an  upper  bound  valid  for  tubes  of  large  radii. 
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